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ABSTRACT 

Gluon distributions in the carbon and tin nuclei are investigated by using a Q 2 rescaling 
model with parton recombination effects. We obtain strong shadowings in the small x 
region due to the recombinations. The ratio Ga{x)/Gn(x) in the medium x region is 
typically 0.9 for medium size nuclei. At large x, the ratio becomes large due to gluon 
fusions from different nucleons. Comparisons with recent New Muon Collaboration data 
for Gsn( x )/Gc(x) indicate that more accurate experimental data are needed for testing 
the model. 
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Recent measurements of gluon distribution ratios Gs n (x)/Gc(x) by New Muon Col- 
laboration (NMC) [1] are the first data which could shed light on gluon distributions in 
nuclei [2]. 

Modifications of the structure function F2(x) in nuclei were discovered by European 
Muon Collaboration (EMC effect) [3]. This effect has been an interesting topic in the sense 
that it may provide an explicit quark signature in nuclear phenomena. On the contrary, 
"gluonic EMC effect" is little known. Gluon distribution functions in the nucleon [4] have 
been investigated by using muon (electron, or neutrino) deep inelastic scattering data [5,6], 
direct photon data [7,8], and muon induced J/ip production data [9,10]. Although there are 
some theoretical predictions [11,12,13,14] for gluon distributions in nuclei, we have little 
experimental data. There are direct photon data for proton reactions with nuclear targets 
[15,16]; however, there is no available large p T data in the WA80 case [15]. Accuracy is 
not good enough for extracting a gluon distribution function of the beryllium (Be) nucleus 
from the E706 data [16]. In any case, we do not expect much modification of the gluon 
distribution in light nuclei such as Be. 

The NMC analyzed inelastic J/ip production data by a color singlet model [17] and 
obtained gluon distribution ratios Gs n (x) /Gc{ x ) [!]• These are interesting data which 
indicate modifications of the gluon distribution in nuclei. We should note that it is very 
important to know the gluon distributions in nuclei. For example, J/ip suppressions in 
heavy-ion collisions were proposed as a signature of quark-gluon plasma [18]. Although 
other initial and final state interactions may explain the J/ifi suppression phenomena, it is 
important that initial conditions as local (gluonic) EMC effects [19] should be subtracted 
out for investigating the physics origin of the suppressions [20] . Such issues have not been 
well studied yet because there is little data for gluon distributions in nuclei. 

As a model for explaining the EMC effect, we take a Q 2 rescaling model [21] and 
apply it for the gluon distribution. The model was proposed originally for the structure 
function i^cc) by considering the possibility that an effective confinement radius for quarks 
is changed in a nuclear environment. In fact, strong nucleon overlaps are expected in nuclei 
by noting the fact that the nucleon diameter is approximately equal to the average nucleon 
separation in nuclei. From this confinement-radius change and the Q 2 evolution equation, 



2 



the Q 2 rescaling model was proposed [21]. In this research, we use the above simple picture 
also for gluons. Namely, the nuclear gluon distribution function in the rescaling model is 
given by Ga{x, Q 2 ) = Gn(x, £, a Q 2 ), where £ A is called as the rescaling parameter. 

Other new features which exist if a nucleon resides in a nucleus are parton recombina- 
tions (fusions) . The mechanism has been investigated for explaining the shadowing region 
F 2 (x < 0.1) in parton models [11,12,22]. The recombinations also have important effects 
on the gluon structure function due to processes shown in Fig. 1. We find that there are 
strong shadowings in the small x region (x < 0.05) and strong anti-shadowings in the 
medium x region [12]. In Ref. 22, it was shown that a model of the Q 2 rescaling with 
the recombination can explain experimental data of F2(x) fairly well in the wide x region 
(0.005 < x < 0.8). 

In this report, the gluonic EMC effect is investigated by the Q 2 rescaling model 
combined with the parton recombinations in Ref. 22. Using this model we calculate gluon 
distributions in the carbon (C) and tin (Sn) nuclei. Results in our model are compared 
with the recent NMC ratios Gs n (x)/Gc(x). 

In order to investigate recombination effects, we calculate contributions to G(x) from 
processes in Fig. 1 as investigated by Close, Qiu, and Roberts [12]. For example, a gluon 
with momentum fraction x is produced by a fusion of gluons from the nucleon 1 and 2. A 
modification of a parton distribution P3(xs), due to the process of producing the parton 
Ps with the momentum fraction x% by a fusion of partons p\ and p2, is given by [12,22] 

Ap 3 (^3) = K J dx 1 dx 2 pi(xi) p 2 {x 2 ) r pi p 2 _ > p s (xi,X2,X3 = Xi + x 2 ) S(x,x 1 ,x 2 ) , (1) 

where K is given by K = 9A 1 / 3 a s /(2i?Q<5 2 ) with i?o=l.l fm and the strong-interaction 
coupling constant is a s (Q 2 ) = 4n/[9 ln(Q 2 / A 2 )]. The 5 function is given by 8{x — x\ — x 2 ) 
for the processes in Figs, la, Id, and le and 8{x — x\) for Figs, lb, lc, If, lg, lh, and 
li. The parton fusion function F(xi, x 2 , X3) is a probability for producing a parton p^ 
with momentum fraction x% by a fusion of partons p\ and p 2 with momenta x\ and x 2 
respectively. 

Now, we discuss numerial analysis. In evaluating recombination contributions AG(x) 
by using Eq. (1) [23], we assume that a leak-out parton (we denote p*{x)) is a sea quark 
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or a gluon and that the momentum cutoff function [12,24] for this parton is taken as 
w(x) = exp(—m 2 N ZQX 2 /2) with zo=0 or 2 fm. Then, distributions for the leak-out partons 
are q*(x) = w(x)q sea (x), q*(x) = w(x)q(x), and G*(x) = w(x)G(x). Input parton 
distributions are given by a recent parametrization in Ref. 4. Q 2 = 5 GeV 2 is used in the 
parametrization and for calculating K in Eq. (1). The QCD scale parameter A in a s (Q 2 ) 
is taken as A=0.2 GeV. In our theoretical analysis, targets nuclei are assumed as 12 C or 
118 Sn. The rescaling parameters for these nuclei are taken from Ref. 21, and they are 
UC)=1.60 and £ A (£ra)=2.24. 

Before discussing the gluon distributions in C and Sn, we first check that our model 
can explain structure functions F^(x) of these nuclei. As shown in Ref. 22, our results are 
consistent with SLAC, EMC, and E665 experimental data for F^(x) of C, Ag, Sn, and Xe 
nuclei in the wide x range (0.005 < x < 0.8). In explaining these data, the most important 
factor is the gluon shadowing. Taking modified parton distributions due to the recombina- 
tions at small Q 2 , we should calculate distributions at large Q 2 , where the structure func- 
tions were measured. Instead of solving the evolution exactly, we simply used a solution 

x d 

[25] for the Altarelli-Parisi equation in the small x region (xSqf ea (x) = — — — [xAG(x)], 
where i=u, d, or s). This approximate way of treating the evolution violates the momen- 
tum conservation even though it is satisfied in the recombinations. The effect due to this 

extra term is given by 6 / dx x Sqf ea (x)= +^ lim x 2 AG{x) + ^ / dxxAG(x). The first 

J 2 x->o 2 J 

term vanishes if the input gluon distribution satisfies, for example, lim xG(x) =constant. 

x— >o 

Dominant contributions to AG(x) come from the gluon-gluon fusion processes; however, 
they satisfy the momentum conservation by themselves (f dxxAGGG->G{ x ) = 0). There- 
fore, the violation is a small effect due to qq — > G, Gq — > q, and Gq — > q processes. A 
numerical evaluation for the Ca nucleus indicates that such violation effect is less than 1% 
(6 / dxx5qf ea {x) = —0.005); hence it is not a serious effect on the momentum conservation. 

Using our model, which can explain at least the structure function F^ix), we predict 
Gc(x) and Gsn(x)- We take the same rescaling parameter for all partons by taking a 
naive consideration that the confinement radius change modifies all parton momenta at 
the same rate. Predicted gluon distributions for C and Sn are shown in Figs. 2a and 2b. 
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In these figures, the dashed curves are recombination effects shown by 



^(recombination) = 1 + -pr^r^T^ > ( 2 ) 



AGa(x, Q 2 ) is calculated by using Eq. (1) (note that no Q 2 rescaling is used) and explicit 
expressions are given in Ref. 22. Solid curves are combined contributions from the rescaling 
and recombinations and they are shown by the ratio 

G A (x, Q 2 ) G A (x,Q 2 ) + AG A (x,Q 2 ) 

G N (x,Q 2 ) ~ G N (x,Q 2 ) • [ } 

In these equations, Ga{x, Q 2 ) and AGa{x, Q 2 ) are given by the rescaling model, Ga{x, Q 2 ) 

=Gn(x, Z, a Q 2 ) and AGa(x, Q 2 ) = AG{x, Z, A Q 2 ). The recombination mechanism produces 
strong shadowing effects in the small x region due to gluon-gluon and gluon-quark fusion 
processes in Fig. 1. In the medium-large x region, contributions are dominated by the 
gluon-gluon fusion process in Fig. la. It is interesting to note in our model that gluon 
distributions at x > 1 could be produced in the fusion process. This is the reason why 
the ratio goes to infinity at x — > 1 in Figs. 2a and 2b. The Q 2 rescaling contributions are 
opposite to the recombination. The rescaling effects are positive in the small x region and 
are negative in the medium-large x region. Combined contributions shown by the solid 
curves in Figs. 2a and 2b indicate strong shadowings in the very small x (x < 0.02) region, 
depletions (0.8—1.0) in the medium x (0.2 < x < 0.6), and large ratios in the large x 
(x > 0.7). 

In investigating shadowings in F2(x) in Ref. 22, we used the rescaling for parton 
distributions at i < 0.1. Although the rescaling produces large positive contributions, 
they are counterbalanced by the large shadowings produced through gluon distributions 
(5F2 in Ref. 22). Therefore, combined contributions are not very dependent (about 5% 
differences) whether or not the rescaling is used in the region (0.005 < x < 0.1). On the 
contrary, gluon distributions at small x are very sensitive to whether or not the rescaling 
is used as shown in Figs. 2a and 2b. If there is no rescaling at x < 0.1, we should have 
strong gluon shadowings in the region x < 0.1 as shown by the dashed curves. However, 
if the rescaling is used, gluon distributions are shadowed only in the very small x region 
(x < 0.02) as shown by the solid curves. 
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The nuclear gluon distributions in the medium-large x region are very sensitive to the 
momentum cutoff as shown in Figs. 2a and 2b. For example, Gs n ( x ) /Gn(x) = 1-02 (for 
zq = 0) at x = 0.4, but it is 0.87 if zq = 2 fm. Because Gn(x) itself is very small at 
x = 0.4, it may seem to be an insignificant problem. However, it is important for describing 
e.g. p T dependence of J/ip in heavy-ion collisions. The rapid increase of Ga(x)/Gn(x) 
in the medium x region could be responsible for the p T dependence of J/ip suppressions 
observed by NA38 [26] , although the p T slope obtained in the local gluonic EMC effect is 
rather small [20] compared with the NA38 data at this stage. The gluon distributions at 
the medium x are so sensitive to the momentum cutoff w(x) that we need to study more 
about the cutoff [24] . We leave the problem of the cutoff as a future research topic. 

Calculated gluon distributions are compared with the NMC data for Gs n ( x )/ Gc( x ) 
[1] in Fig. 3. The dashed curve shows recombination results with zo=0 and the solid 
(dash-dot) curve shows combined results of the Q 2 rescaling and the recombinations with 
zo=0 (zq = 2 fm). Our theoretical results shown by the solid and dash-dot curves indicate 
ratios Gs n (x)/Gc{x) < 1 at x < 0.03 due to the gluon shadowing. The ratios are 
about 1.03 in the region (0.05 < x < 0.20) and are dependent on the momentum cutoff 
w(x) in the medium x (x > 0.2). We notice that experimental errors are very large in 
comparison with typical theoretical modifications. Because our model predictions for the 
modifications are less than 5% in the range 0.05 < x < 0.2, we need accurate data better 
than 5% accuracy in order to test the model. Because the NMC data for Gs n (x)/Gc(x) 
are not accurate enough, we should wait for better measurements of Ga(x), for example 
a proposed experiment at RHIC [27], in the small x region for investigating details of 
the gluon shadowing. It is also important to know gluon distributions in the medium x 
region (0.3 < x < 0.6), although distributions themselves are very small, for studying 
J/ip productions in heavy-ion collisions [20]. 

In summary, we investigated gluon distributions in the nuclei C and Sn by the rescaling 
model with parton recombinations effects. We obtained shadowings in the nuclear gluon 
distributions in the small x region due to the recombinations and depletions [typically 
Ga(x)/Gn(x) ~ 0.9] in the medium x region. The ratio Ga(x)/Gn(x) becomes large 
at x > 0.6 due to gluon fusions from different nucleons. The ratio in the medium-large 
x region is very sensitive to the momentum cutoff for leak-out partons in our model. 
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Comparisons with the NMC data indicate that more accurate experimental data are needed 
for testing the model. 
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Figure Captions 



1. Modifications of G(x) due to parton recombinations. 

2. Predicted gluon distributions in the nuclei (a) C and (b) Sn. The dashed curves 
indicate parton recombination effects (Eq. (2)) and the solid curves are combined 
results (Eq. (3)) of the recombination and the Q 2 rescaling. £ A =1.60 for C and 2.24 
for Sn. 

3. Comparisons of our theoretical results with NMC data [1]. The dashed curve 
indicate parton recombination effects with zq=0 and the solid (dash-dot) curve shows 
combined results of the recombination and the Q 2 rescaling with zq=0 (zq=2 fm). 
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